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L'invention concerne un procede pour estimer la reponse 
impulsionnelle d'un canal de propagation, notamment ses parametres tels 
que les valeurs de retard, les directions d'arrivee, ainsi que les valeurs des 
5 amplitudes complexes associees k ces parametres, avec une connaissance 
a priori du signal. 

Elle s'applique, par exemple, pour estimer les parametres d'un 
filtre a reponse impulsionnelle finie, pouvant mais non hecessairement, 
s'ecrire de fagon speculaire, c'est-a-dire que le filtre peut s'ecrire sous la 
10 forme de Dirac ponderes. 

L'invention est utilisee pour estimer des canaux de propagation 
dans le domaine des radiocommunications, mais aussi de maniere generate, 
a tout signal filtre par un filtre lineaire a reponse impulsionnelle finie. 

15 Dans un systeme de transmission, notamment par ondes radio, un 

emetteur emet un signal dans un canal, de transmission a destination d'un 
recepteur. Le signal emis subit des fluctuations d'amplitude et de phase dans 
le canal de transmission, le signal regu par le recepteur est constitue de 
copies temporellement decalees et modifiees du signal emis. Les fluctuations 

20 du signal et les decalages generent ce que THomme du metier appejle 
('interference intersymbole. Linterference provient notamment de la loi de 
modulation employee pour la transmission et egalement de la propagation de 
multi-trajets dans le canal. 

Le signal re?u est generalement issu d'un grand nombre de 

25 reflexions dans le canal, les differents trajets empruntes par le signal emis 
conduisant ainsi a des retards varies au niveau du recepteur. La reponse 
impulsionnelle du canal represente Pensemble des fluctuations auxquelles 
est soumis le signal emis. 

L'estimation du canal de propagation dans un systeme de radio- 

30 communication est interessant a plusieurs egards, dont certains sont 
indiques ci-apres a titre d'exemples. 

• Les demodulateurs necessitent generalement la connaissance des 
canaux afin de remedier aux effets nefastes qu'ils ont induits, 
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• Un deuxieme interet est la localisation urbaine ou extra-urbaine par 
exemple le principe de localisation d'urgence « 91 1 » aux Etats-Unis, 

• Enfin, la connaissance des canaux de propagation peut egalement servir 
a ['utilisation d'antennes intelligentes (de type smart antennas) tant a la 

5 reception qu'a remission. 

Differentes techniques sont connues de Tart anterieur pour estimer 
les canaux de propagation et leurs parametres. 

Par exemple, le document de R. Rick and L. Milsteil intitule 
"Performance acquisition in mobile ds-cdma systems", IEEE Trans on 
10 Communications, Vol: 45(No: 11):pp: 1466-1476, November 1997, propose 
de rechercher les retards de propagation en utilisant un banc de detecteurs 
non-coherents. Les resultats sont proposes pour des canaux a multi-trajets, 
en presence de Doppler et d'interferences inter-cellules et intra-cellulaires. 

Le document de R. Rick et L. Milsteil ayant pour titre "Optimal 
15 decision strategies for acquisition of spread spectrum signals in frequency 
selective fading channels" paru dans IEEE Trans, on Communications, Vol: 
46(No: 5):pp: 686-694, May 1998, divulgue une regie de decision optimale 
en se fondant sur les sorties des correlateurs proposes dans le document 
precedemment cite. Une telle technique mono-utilisateur est limitee par les 
20 interferences dans le cas de multi-utilisateurs. 

II est aussi connu d'utiliser des filtres de mise en forme 
rectangulaires, par exemple selon des precedes decrits dans Tun des 
documents suivant: 

• E. Strom, S. Parkvall, S. Miller, and B. Ottersen intitule "Propagation 
25 delay estimation in asynchronous direct-sequence code-division multiple 

access systems", IEEE Trans on Communications, Vol: 44:pp: 84-93, 
January 1996 

• S. Parkvall intitule "Near-Far Resistant DS-CDMA Systems : Parameter 
estimation and Data Detection", PhD thesis, Royal Institute of Technology 

30 Stockholm, Sweden, 1996. 

• S. E. Bensley and B. Aazhang intitule "Maximum likelihood estimation of a 
single user's delay for code division multiple acess communication 
systems", Conf. Information Sciences ans Systems, 1994. 
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Dans le cas de filtres de mise en forme de duree superieure a un 

♦ 

temps de chip, ces differents procedes ne conviennent plus. 

Des algorithmes d'estimation conjointes des angles d'arrivee et 
des temps de retards differentiels, sur des signaux recus et connus, fondes 
5 sur des techniques de sous-espaces ont ete proposees par exemple dans le 
document de P. Gounon intitule "Analyse spatio-temporelle haute resolution 
a I'aide d'une antenne active", Traitement du Signal, Vol: 1 1 (No: 5):pp: 351- 
.360, 1 994. 

Le document de A. J Van der Veen, M. C. Vanderveen, et A. J. 

10 Paulraj, ayant pour titre "Joint angle and delay estimation using shift- 
invariance properties", IEEE Sig. Proc Letters, Vol.4(No.5):pp:142-145, 
1997, divulgue des procedes d'estimation des parametres physiques de la 
propagation a I'aide de methodes a base de sous-espaces. 

De telles methodes souffrent toutefois d'une degradation des 

15 performances des que les reponses impulsionnelles des canaux de 
propagation sont correlees. Cette situation se produit notamment lorsque ies 
amplitudes complexes ne varient pas assez rapidement sur la matrice de 
covariance des reponses impulsionnelles estimees au sens des moindres 
carres a I'aide du signal emis. 

20 La figure 1 montre differentes techniques du Maximum de 

vraisemblance pouvant etre utilisees en fonction des hypotheses. 

Une methode au sens du maximum de vraisemblance a par 
exemple ete proposee dans Tune des references suivantes : 

• J. Grouffaud, "Identification spatio-temporelle de canaux de propagation 
25 a trajets multiples", PhD thesis, Ecole Normale Superieure de Cachan, 

Juin 1997. 

• M. Wax and A. Leshem, "Joint estimation of delays and directions of 
arrival of multiple reflections of a known signal.", IEEE Trans, on Signal 
Processing, Vol: 45(No: 10):pp: 2477-248, October 1997. 

30 mais ne traitent pas le contexte MIMO (abrege anglo-saxon de Multiple Input 
Multiple Output). 
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Le document de P. Graffouliere intitule "Methodes actives spatio- 
temporelles large bande", publie dans Techniques et performances. 
Applications en sonar. PhD thesis, I'lNPG, March 1997,divulgue egalement 
['estimation au sens du maximum de vraisemblance mais les etudes de 
5 performances ne traitent que le cas d'une source unique ou de plusieurs 
sources tres nettement separ§es. Une etude similaire est divulguee dans le 
document de N. Bertaux ayant pour tit re "Contribution a I'utilisation des 
methodes du Maximum de Vraisemblance en traitement radar acti", PhD 
thesis, Ecole Normale Superieure de Cachan, Janvier 2000, Pour des 

10 applications de radars actifs dans le cas monosource. 

La presente invention concerne un procede qui permet, 
notamment, d'estimer les parametres d'un canal de propagation, en 
travaillant sur les signaux correles, en selectionnant un certain nombre 
d'echantillons et en recherchant les valeurs des parametres retard et/ou 

15 directions d'arrivees, par exemple, qui permettent de reconstruire au mieux le 
signal regu. 

Elle concerne egalement un procede integrant des techniques de 
compression d'impulsion comme pretraitement. 

20 L'invention concerne un procede pour estimer un ou plusieurs 

parametres d'un canal de propagation avec une connaissance a priori du 
signal dans un systeme comportant un ou plusieurs capteurs. 

II est caracterise en ce qu'il comporte au moins les etapes 

suivantes : 

25 • correler le signal regu par le ou les capteurs s(t) avec un signal connu 
c(t), 

• selectionner un nombre d'echantillons du signal obtenu apres I'etape de 
correlation, 

• determiner les parametres du canal tel que x et/ou 0 qui permettent de 
30 reconstruire au mieux les signaux regus en utilisant une methode au 

maximum de vraisemblance. 
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Selon un mode de realisation, les caracteristiques du systeme de capteurs 
sont connues et il comporte par exemple : 

• une etape de correlation avec un signal connu c(t) egal a 1 , 

• les signaux recus sur I'antenne s'exprime sous la forme X=S(9,-c)h+B 
5 • les estimees des parametres x et 6 s'expriment sous la forme : 

1 * 
i _ . ♦ .. . * ' ' m 

• • v- • 

6,t = argroin||n£(0,-r)X|| 2 
= argmin{X*n^(0 > T)X} 

0 • ob n x s est le projecteur orthogonal a I'image engendree par les vecteurs 

colonnes de S(9,i:). 

Les caracteristiques du systeme de capteurs sont par exemple la 

reponse a I'antenne. 

II peut comporter une etape de determination des amplitudes 
1 5 complexes h de la reponse impulsionnelle du canal de propagation a partir 
des estimees des parametres estimes % et 9. 

Selon un autre mode de realisation les caracteristiques du systeme de 
capteurs ne sont pas connues et il comporte par exemple : 

• une etape de correlation des signaux re?us par le reseau de capteurs 
20 avec un signal connu c(t) egal a 1 , 

• une etape ou les signaux re?us s'expriment sous une forme concatenee 
Y =*F(i;)a+N ou ^(t) est egal au produit convolue de la matrice unite l N 
avec la matrice S{r) = [s'tf), -- .s 1 ^), •• ,s(t^)] et a contient les 
reponses des trajets des differents utilisateurs, 

25 • une etape d'estimation du vecteur retard i a partir de 

t = argmin ||ni(T)Y|| 2 

T 

= argmintr(Y t ni(r)Y) 



1er depot 




6 



ou n x * est le projecteur orthogonal a Pimage engendree par les vecteurs 
lignes de 

L'etape de correlation des signaux est effectuee avec un signal 
c(t) different de 1, en ce que les caracteristiques du systeme de capteurs 
5 sont connues et en ce quMI comporte une etape d'estimation des parametres 
x et 6 a partir de 

> Ui=l- r) (*♦(*, t)!^ 1 *^, r))" 1 

■ • - 

10 

LMnventipn conceme aussi un dispositif pour estimer un ou 
plusieurs parametres d'un canal de propagation avec une connaissance a 
priori du signal dans un systeme comportant un ou plusieurs capteurs. 

II est caracterise en ce qu'il comporte au moins : 
15 • un dispositif adapte a correler le signal re$u par le ou les capteurs s(t) 
avec un signal connu c(t), 

• un dispositif adapte a selectionner uh nombre d'echantillons du signal 
obtenu apres l'etape de correlation, et 

• uh dispositif adapte a determiner les parametres du canal par une 
20 methode au maximum de vraisemblance. 

Le procede est applique par exemple dans des applications de 
type MIMO (Multiple Input Multiple output) ou SIMO (Single Input Single 
Output). 

25 

(-'invention presente notamment comme avantages : 

• dans le contexte MIMO cooperatif, Festimateur utilise dans le cadre d'une 
antenne connue, utilise un modele parametrique des signaux re?us sur 
I'antenne, contrairement au maximum de vraisemblance deterministe qui 

30 estime chaque echantillon du signal sans a priori, 
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• en lirnitant le domaine d'analyse dans certains cas d'application, elle 
conduit a une reduction de la taille des donnees traitees d'ou un 
traiternent plus rapide. 

D'autres caracteristiques de Tinvention apparaitront mieux a la 
lecture de la description qui suit ou sont proposes des exemples de mise en 
oeuvre a titre illustratif, ceci en reference aux figures annexees qui 
represented : 

• La figure 1 represente differentes approches d'estimation parametriques 
au sens du Maximum de vraisemblance, 

• La figure 2 schematise une premiere variante du procede selon I'invention 
pour laquelle les caracteristiques des antennes de reception sont 
connues, 

• La figure 3 represente un schema d'une deuxieme variante de realisation 
ou les caracteristiques des antennes de reception ne sont pas connues/ : 

• La figure 4 schematise une autre variante du procede comportant une 
etape prealable pouvant s'appliquer aux methodes decrites aux figures 2 
et 3 precedemment. 

• 

Afin de mieux faire comprendre I'objet de I'invention, la description 
est donnee a titre illustratif et nullement iimitatif dans le domaines des 
radiocommunications mobiles, dans un contexte cooperatif ou le recepteur 
dispose de certaines informations, par exemple, la connaissance dans une 
certaine mesure du signal emis (connaissance a priori du signal). 

Le recepteur est par exemple constitue de N capteurs identiques 
et omnidirectionnels. Les signaux s(t) re9us par chacun des capteurs etant 
echantillonnes regulierement a la cadence T e . Cet exemple, n'empeche pas 
1'application du procede selon Finvention a d*autres types de recepteurs. 

Avant de detailler les differentes variantes de mise en ceuvre du 
procede, quelques rappels sur la moderation des signaux et quelques 
hypotheses permettant la mise en ceuvre du procede sont enonces. 

Modelisation des signaux regus sur le systeme de reception 

En presence de U utilisateurs (sources), les U signaux s u (t) avec 
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1 < u < U se propagent a travers differents canaux avant d'etre re9us par un 
ensemble de capteurs indices k. La reponse impulsionnelle du canal entre 
I'utilisateur u et le capteur k est notee h u .(t,t), ou t est I'instant d'arrivee du 
signal et t son retard. 

En appelant h u k (t,t) le vecteur de dimension {Nxl) dont les 
composantes sorit les h u k (t,T), avec 1 < k < Ket K le nombre de capteurs, on 
peut ecrire le modele des signaux re$us sous la forme (1) : 

x(t)=2s B (T)*h u (t,t)+b(t)ou 
b(t) est le vecteur bruit cie dimension (Nxl). 

Hypothese speculaire 

Temporellement speculaire 

Le modele speculaire est un modele de canaux largement utilise 
pour lequel, on considere que les trajets sont vus par Pantenne comme etant 
des trajets temporellement ponctuels. 

La fonction de transfert du canal devient done (2): 

h u (t,x)=2 h u p (TpU(t))5(t.T p u(t)) 

avec 

• p u le nombre de trajets de propagation de Tutilisateur u , 

u 

• n p (T p u (t)) le (Nxl) vecteur contenant le produit des reponses de 
I'antenne au trajet p par I'amplitude complexe du trajet. 

Le signal regu s'ecrit dans le cas du canal speculaire (3) 

x (*) = E 5Z h P - r p u (t)) + h(t) 
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Antenne calibree 

Si I'antenne est calibree, avec I'hypothese que les signaux arrivant 
sur I'antenne ont une dispersion angulaire negligeable et que I'hypothese 
bande etroite des signaux est respectee, le signal recu est modelise par (4) : 



x to - s £ - t p + b(<) 

u=l p=l 

oD 

10 • e p u est la direction d'arrivee du trajet p de I'utilisateur u, 

• h " n est le gain complexe du trajet p de I'utilisateur u recu sur le 

premier capteur. 

L'hypothese precedemment faite (specularite des directions 
d'arrivees) n'est pas restrictive dans le cas d'une station de base en hauteur. 
15 si la specularite n'est pas suffisante, le modele peut etre modifie en ajoutant 
une largeur de lobe de direction d'arrivee qui sera estime par la methode 
proposee selon I'invention en ajoutant les parametres de ce lobe dans les 
parametres a estimer. 

20 Les differentes variantes du procede selon I'invention, font appel a 

la methode du maximum de vraisemblance, dont differentes approches 
d'estimation parametrique ont ete donn6es a la figure 1 . 



Premiere variante de realisation - antenne calibree - 

25 Une premiere variante de realisation du procede selon I'invention 

utilise une methode du maximum de vraisemblance adaptee aux signaux 
emis, dans le cas ou les signaux emis sont en partie connue, et la reponse 
de I'antenne pourtoutes les directions d'arrivee est connue. 

La figure 2 represente un schema bloc ou le bloc de traitement t , 

30 comportant par exemple un microprocesseur, adapte a mettre en osuvre le 
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procede selon I'invention, recoit d'une part les differents signaux recus par 

les antennes 2i, ainsi qu'un signal connu correspondant au signal emis. 

En sortie du bloc de traitement le procede selon I'invention 

dispose au moins des estimees des directions d'arrivee G et des vecteurs de 
5 retard tqui sont par exemple utilisees pour en deduire la reponse 

impulsionnelle du canal h(t), ceci pour une utilisation donnee ou pour 

I'ensemble des utilisateurs. 

Les signaux re9us sur le reseau d'antennes sont des versions 

retardees, dilatees et eventuellement transposes en frequence des signaux 
10 emis par les utilisateurs presents dans une cellule (ou un canal de 

propagation). Ces signaux s'expriment done comme une combinaison 

lineaire d'un ensemble de fonctions de base connues et dont on veut extraire 

les parametres. 

Le procede est presente dans le contexte general d'estimation 
15 parametrique des angles 0 et des retards i pour une antenne quelconque 
constitute de N capteurs identlques et omnidirectionnels. Les signaux sont 
echantillonnes regulierement a la cadence T e . On suppose le nombre de 
trajets connu et, si tel n'est pas le cas, on le determine en utilisant une 
methode de detection du nombre de sources qui permettent d'obtenir le 
20 nombre de trajets, par exemple decrites dans I'une des references suivantes: 

• "Estimation of the number of signals in the presence of unknown 
correlated sensor noise" de J J.Fuchs, IEEE Trans.on signal Processing, 
Vol 40(N° : 5) pages 351-360, 1994; 

• "Methodes a Haute resolution pour le traitement d'antenne et I'analyse 
25 spectrale" de J.Gouffraud, P.Larzabal et H.CIergeot, chapitre 6, Selection 

d'ordre de modele, aux Edistions Hermes, Fevrier 1 998; 

• "A new look at the statistical model identification" de H.Akaike, IEEE 
Trans. On Automatic Control, Vol 19(N°6), december 1974, 

• "Detection of signals by information theoric criteria" de M.Wax et 
30 T.Kailath, IEEE Trans.on ASSP, Vol 33(N°2), pages 387-392, April 1985. 



20 
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Les etapes mises en ceuvre dans I'algorithme selon I'invention 
restent valables lorsque le nombre de sources est surestime. 
Modelisation 

Le signal x(t) (equation (4)) recu par le reseau de capteiirs 2i est 
5 echantilione au ryhme 7 e> pendant une duree f=[Te,A/ e r e ] avec Ne nombre 
d'echantillons ou nombre d'observations pendant laquelle les amplitudes 
complexes de la reponse impulslonnelle pour un utilisateur u et un trajet p, 
h tt (tp u (t)) sont consideres comme constantes. 

Ces N e observations de dimension Mx1 sont concatenees dans 
10 le vecteur X de dimension A/©Ak1 : X=[ x(T e ) T , -,x(N e T e ) T ] T (5) ou 
I'exposant T designe la transposee. 

On definit le vecteur B de dimension WgNxl contenant la 
concatenation des echantillons du bruit durant la meme periode par (6) : 

15 B=[b(T e ) T ,",b(N e Te) T ] T 

Enfin nous notons s u (t p u ) le vecteur de dimension A/gxl 
correspondant au signal emis par I'utilisateur u, retarde de t p u et 
echantillonne : 



s u (V)4 s U <Je-0. -.s u (N e Te- x P u )f (7) 



25 



Le vecteur des concatenations des signaux re?us peut done s'ecrire (8) 



x = EEWW) ® < e W + B 



ou ® est le produit de Kronecker. 
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V 

Notons S(9,t) la mat rice N e NxP ( P=X pU ) des vecteurs 
retard-direction definie de la maniere suivante (9) : 

5 S(0,t)= [[s 1 ^ 

OU 

- 9 est le vecteur reel de taille P des directions d'arrivee 
10 - x est le vecteur reel de taille P des retards, (P etant le nombre de 
trajets) 

On peut done ecrire Tequation des signaux re9us sur Pantenne 
comme : X=S(9,x)h+B (10) ou h est le Px1 vecteur complexe des 
amplitudes. 

15 Cette formulation prdsente I'interet d'exprimer les observations 

comme un modele lineaire fonction des amplitudes complexes h et du bruit. 

Notons bien qu'il s'agit ici d'une moderation MIMO (Multiple input 
multiple output), regroupant I'ensemble des signaux de tous les utilisateurs. 

II n'est pas necessaire de connaTtre tout le signal emis. par 
20 exemple en communications mobiles, sequence de symboles connus sont 
emis periodiquement resultant en un signal connu. En radar ou en radio- 
altimetrie, par exemple, I'impulsion emise est connue. 

Log-vraisemblance 

Les amplitudes complexes h sont considerees comme inconnues 
25 mais deterministes. Dans le but de simplifier Pexpose de invention sans 
toutefois la limiter, les signaux s u (f) formant S(9/t) sont consideres comme 
parfaitement connus aux parametres 9,x pres. Cela conduit a un modele 
d'observation (1 0 ) ou seul le vecteur bruit additif est a!6atoire, de densite de 
probability gaussienne, represents par (11) : 
30 / 

P{X|<r 2 ,h,V} = nN J a2NeN exp (-^ (X - 8(0, r)h)> (X - S(0, r)h) j 
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La log-vraisemblance des observations est dont (12) : 



1 



L (X|<7 2 , h, 6, T) = -N e N l0g(7T<T 2 ) - ^ ||X - S(0, r)h|| 



Estimateur des parametres <y*,Y\,9,% 

Grace a la separabilite des parametres o 2 , h, 9, x, I'expression de la 
log-vraisemblance (12) autorise la determination analytique du vecteur des 
10 amplitudes complexes h et de la puissance du bruit o 2 en fonction de 9 et x 
Ces estimees sont toutes deux obtenues en annulant les derivees de la log- 
vraisemblance respectivement par rapport a o 2 et par rapport a h . Elles ont 
pour expression respectivement (13) et (14): 



15 



a 2 = 



HX-S^t)!^ 



et 



h= {sHe^sie^))- 1 s\e,r)x 



20 



En substituant les valeurs de h et de a 2 dans la log-vraisemblance 
(12), I'expression se simplifie et les estimees de 9 et x sont donnees par (15) 



25 



0,t = argmin 

0,T 



x - s(0,r) {s\e,T)s(e,r))- 1 s'(0,t)x 



30 



En notant n£ le projecteur orthogonal a I'image engendree par les 

* 

vecteurs colonnes de S(9, x) (16) 

U £(0, t) = I — 8(0, t) (S* (0, r)S{e, r)) _i S* (6 , r) 



1 er depot 





14 



Les estimees sont les solutions de (17) 

■ 

■ m * 

0,9 = argmml|ni(0,T)X|| 2 
5 = argminlXtni^.rJX} 

Ces vaieurs estimees sont ensuite utilisees respectivement dans 
les expressions (13) et (14) pour determiner la puissance du bruit et la 
reponse impulsionnelle du canal de propagation. 
10 Le critere introduit dans (1 7) est note dans la suite du texte par 

L(0 t r) = Xtn£(^T)X 

Mise en oeuvre du procede selon l'invention dans le cas d'une 

15 antenne non calibree 

La figure 3 schematise un exemple de mise en oeuvre du procede 
selon invention dans le cas ou les caracteristiques du recepteur ou de 
Tantenne ne sont pas connues, Tantenne est formee de plusieurs capteurs, 
dont la reponse n'est pas connue designee souvent par I'expression 
20 "antennes non structurees". 

Les anfennes 2i regoivent les signaux qui sont transmis a un bloc 
de traitement 4 adapte a mettre en oeuvre les etapes du procede, et qui 
re9oit lui-meme un signal connu. 

Le bloc de traitement equipe par exemple d'un microprocesseur 
25 programme en consequence delivre au moins la valeur estimee du vecteur 
de retards % et la valeur estimee de la reponse impulsionnelle du canal de 
propagation. 
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Modelisation 

Le signal recu s'exprime sous la forme de I'equation (4) et est 
echantiilonne au rythme Te sur une periode de temps t = [Te, NeTe] dans 
laquelle les amplitudes complexes de la reponse impulsionnelle du canal de 
propagation hp o (x p u (0) peuvent etre considerees comme constantes ou 
sensiblement constantes. Ces Ne echantillons de dimension Nx1 sont 
concatenes dans un vecteur Y de dimension NeNxl tels que (18) : 



Y = 



xi(T e ),--- ,*i(AT ? T«),. • ,x N (T e ),- - ,x N (N e T e ) 



capteur 1 



capteur N 



Ou xi(nTe) est I'echantillon n du capteur i. 

De meme, le vecteur N de dimension NeNxl contient la 

15 concatenation des echantillons de bruit (19) 



N = 



1 T 



bi (T.), , • , , h(NeT c ), • • • , b N (T e ), : . • , b N (N e T e ) 



capteur 1 



capteur N 



20 Ou bi(nTe) est I'echantillon de bruit n du cateur d'indice i. 



Le signal Y peut s'ecrire (20) 



Y=¥0c)<x+N 
25 ou la matrice (N e NxNP) 
comme suit (21) : 



*F est construite a partir des differents s^V) 



30 



ou S(t ) est la matrice (NgxP) telle que (22) : 
Sir) - [s 1 ^ 1 ), • - • • • • , s (r£,)] 
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a contient les reponses des capteurs aux differents trajets des 
differents utilisateurs (23) : 



10 



oc = 



T T 



avec 



capteur 1 



15 Les amplitudes complexes sont toujours supposees constantes 

pendant la phase d'estimation des parametres. 

La moderation des signaux refus donnee par Pequation (20) 
nous fournit un modele affine. Comme precedemment, cette expression 
prend en compte les signaux de Pensemble des utilisateurs. 

20 Log-vraisemblance 

Les amplitudes complexes a sont considerees dans ce cas 
comme etant inconnues mais deterministe. Les signaux qui interviennent 
dans ¥(t) sont supposes connus aux variables x pres. Cela conduit a un 
modele ou seul le bruit N est aleatoire, de composantes gaussiennes et de 
25 densite de probability (25) : 

P {Y\a\ a, t) = ^ N J a2NeN exp (~±tr ((Y - *(r)a)t(Y - *(r)a))) 

30 La log-vraisemblance est alors donnee par (26) 
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L(Y|aV,T) = -^ 

5 Estimateur des parametres o^at 

Comme dans le cas du maximum de vraisemblance adapte aux 
signaux recus avec une antenna connue tel qu'expose precedemment, les 
expressions analytiques des amplitudes complexes a , ainsi que la 
puissance du bruit o 2 peuvent etre exprimees en fonction du jeu de 
1 o parametres des vecteurs retards x. 

Ceci s'obtient par exemple en annulant respectivement les 
derivees de la log-vraisemblance par rapport a a et par rapport ao 2 . 

On obtient I'expression (27) 



15 N R N 



o 2 = ™|jY-^(r)a|| 2 



Et(28) 

a 



= («nr)*(r))"' *V)Y = ¥*(<r)Y 



20 En remplagant a et o 2 par leurs estimees, la log-vraisemblance se simplifie 
et I'estimateur de x est (29) 



t — arg mm 

T 



Y - *(t) (**(t 1 **(t)Y 



25 En notant n A (T) le projecteur orthogonal a I'image engendree par les 



v 

vecteurs lignes de ¥(t) (30): . 



ni(r) = I - *(t) {t))~ 1 &{t) 



30 I'estimee de x est donnee par (31 ) 
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I r = argimn||ni(r)Y|| 2 

= argmintr (Y f ni(r)Y) 

5 

Les deux variantes de realisation qui ont ete precedemment 
exposees traitent des signaux echantillonnes en sortie de capteur ou encore 
pour un signal correle avec un signal cbnnu de valeur 1. Les matrices a 
traiter sont de tailte NeNxP. 
10 Selon une autre variante de mise en ceuvre du procede qui peut 

s'appliquer aussi bien dans le cas d'une antenne calibree que dans le cas 
d'une antenne non calibree, le procede selon Pinvention utilise des signaux 
correles. 

La figure 4 represente un schema bloc dans lequel un correlateur 
15 10 regoit d'une part les signaux regus par les antennes et d'autre part le 
signal connu. 

Le correlateur 10 comporte un nombre de sorties identiques au 
nombre d'entrees, chaque signal si(t) arrivant sur une antenne etant correle 
par un signal connu c(t). Les differents signaux correles si(t)*ci(t) ou i est 

20 Tindice du capteur, sont ensuite transmis a un dispositif 11 adapte a 
selectionner un nombre d'echantillons, afin notamment de travailler sur un 
intervalle de temps prealablement defini- Les signaux correles et appartenant 
a Tintervalle de temps retenu, sont ensuite transmis a un bloc de traitement 
12 comportant par exemple un microprocesseur adapte a estimer les valeurs 

25 9 et x, lorsque les caracteristiques de Tantenne sont connues ou la valeur de 
x pour une antenne ayant des caracteristiques inconnues. A partir de ces 
valeurs il est possible de determiner Festimee de la reponse impulsionnelle 
du canal de propagation. 

Le correlateur 10 peut etre de type analogique ou numerique, 

30 Tetape d'echantillonnage etant reaiisee avant ou apres Tetape de correlation. 
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Le bloc de trartement regoit aussi le signal emis connu lui-meme 
correle par exemple au moyen d'un correlateur 13. 

Lorsque Ton utilise un intervalle de temps connu, il est utile de 
connaTtre Tinstant de reception du bloc ainsi que la longueur du canal. 

Les signaux consid6res au niveau du precede correspondent aux 
signaux en sortie de la correlation par le filtre adapte au signal emis. 

Procede pour une antenne dont les caracteristiques sont 

connues 

Dans cette partie on utilise les memes hypotheses que dans la 
partie precedente exposee dans le cas d'une antenne dont on connait les 
parametres. On considerera toujours les signaux comme connus et on ne 
s'interessera qu'a Tutilisateur a> pour simplifier la presentation. Ce modele. 
se generalise immediatement au contexte MIMO (Multiple input Multiple 

* 

Output). 



Modelisation 

En correlant le signal regu sur chaque capteur avec le signal.de 
I'utilisateur co , on obtient (32) : 

yw(r > = {^Si*^^ 3 ^ ~ & + b(<) ) * -*■(-*) 

jpll * 



En limitant les donnees traitees a la partie "contenant" le canal de 
propagation en sortie de la correlation, on reduit grandement la complexity 
car on diminue la taille des matrices utilisees. 
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En concatenant les echantillons fitres en sortie des capteurs sur 
un intervalle de temps Nl, le signal X peut etre modelise de la maniere 
suivante (33) : 

X = #(0,<r)h + b 
X= [^'(T.),-.. ,y 0, (JV,r,)] T 



10 



ou 



et 



avec 
et 



h w = 



{(h(t)+b'{t)) T *s»*(-t)} 



$ p =[$J(T e ),-,$^T e )] T 
0) = (s"(t - r p u ) * s"* (-*)) ® a(0) 



15 Le bruit b* 0 suit une loi normale de covariance a?Rb ou la matrice 

R b =R s M i N . R s ° est la matrice de covariance du vecteur (34): 



20 



s w = f s -(r e ),.,. ,5»(Ai7;)jt 



La vraisemblance du signal X* est donnee par (35) 



25 



P(X"|*(*;x),h«->§.,R b ) = 



L___ pvn (-i( X "-*( e ^)h W ) , Rb'(X---*(0,x)h")) 

7T^4^|R b | eXpV ' ' 



d'ou la log-vraisemblance (36) : 



30 



L(X w |*(0,T) > h w ,<4,R b ) = -M^log(7raf.)-log(|R b |)— 1 jlR^ (X" - $(0,T)h" 



En annulant les derivees de la log-vraisemblance par rapport a o 2 
et par rapport a h* 0 , on obtient les estimateurs de et de (37) : 
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a T = 



R b s (X w -*(0iT)h w ) 



NiN • 



et 



h? = ($t(0,T)*(0,T)) r)XT 



En rempla9ant o 2 et h 10 par leur expression dans (36 ) on obtient I'estimateur 
recherche pour les parametres spatio-temporels 8 et % (38) ; (39) 

0, t = arg miri X"R^Ili(0, t)X* 

OU 



Les paramdtres spatio temporels estimes sont par exemple 
utilises pour determiner I'estimee de la reponse impulsionnelle du canal de 
propagation. 

Antenne non calibree 

Dans cette partie, nous nous int6ressons a une methode du MVA 
reduit avec une antenne inconnue. La demarche utilisee sera ideritique a la 
precedente. Nous nous contenterons de donner le modele du signal, la log- 
vraisemblance et les estimateurs. 
Modelisation 

Le signal regu est filtre par sf*(rf) (40) : 
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x"(r) = 



= ( E E - + m*)) * ^ (-*) 



= E h * W ( s "(* - t p ) * ^M)) + (b(t) + b'(t)) * s"(-t) 



5 En concatenant les echantillons en sortie des capteurs puis les 

vecteurs obtenus on obtient (41 ) : 



10 



15 



ou 



et 

avec 
et 



\{(bi(t) +t(t))*s»(-t)\ ) } 

* 

s„ = [s p (2;),... ^(NtTjf 



Le bruit n w suit une loi normale de covariance crlR„ ou la matrice 

i ft 



20 R n =I N ®R sa . 



La vraisemblance du signal Y* est donnee par (42) : 



P(Y"|S(t), <*<*>?., R n ) = 



7r"«"4"'" |R n | 



exp 



^--"-(V -S(r)a" ) 1 (X" -S(x)ot w 



25 



d'ou la log-vraisemblance (43) : 



LCY^r), a", 4,R„) = -JV/Wlogfraf.) - log(|R„|) - -5- 



| R ;5 ( Y "'-§(r)a«*')| 



5 



20 
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En annulant la derivee de la log-vraisemblance par rapport a cr* et par 
rapport a a w . on obtient les estimateurs de cr* (43) et de a M (44) : 



et °r = 



Rn s (Y w - S(r)a w ) 



a" = (SCr)!^ 1 ^))- 1 S(r)R^X 



10 En remplacant dans (42) et par leur expression, on trouve 

un estimateur pour le parametre x (45): 



r = aigminY w R- 1 n^(T)Y w 
ou T 
15 n£ = I - S(t) (S^R-^r)) 'SWR; 1 



La valeur estimee du retard x est ensuite utilises par exemlpe pour 
estimer la valeur de la r6ponse impulslonnelle du canal de propagation. 



Les differentes variantes de mise en oeuvre du precede selon 
I'invsntion, s'applique dans de nombreux domaines du traitement du signal. 
Dans les systemes de communications, il est souvent necessaire de 
connaitre le canal de propagation afin de demoduler les signaux emis. 

25 Compte tenu de la congestion spectrale et de la demande 

croissante des utilisateurs, les systemes future feront communiquer plusieurs 
utilisateurs simultanement dans la meme bande de frequence. De ce fait, les 
techniques proposees ici dans le cadre MIMO ont une application directe. 

De la. meme facon, la technique d'estimation du filtre avec 

30 reduction de la complexity s'applique dans le contexts MIMO et SIMO (Sigle 
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Input Multiple Output). Nous ne revendiquons pas la technique d'estimation 
SISO (Sigle Input Single Output) qui est bien connue. 
Dans le cas d'un traitement radar, ces techniques s'appliquent egalement. 
Elles peuvent etre utilisees dans un contexts multi-cible ou mono-cible avec 
5 echo. 

De meme dans le cadre de la radio-altimetrie, plusieurs echos du 
signal emis peuvent etre requs et tous necessitent d'etre estimes. 

Sans sortir du cadre de Tinvention, le procede peut tenir compte 
10 de Teffet Doppler. Pour cela, une maniere consiste par exemple a estimer les 
vitesses en minimlsant -une fonction multidimensionnelle integrant tous les 
parametres du modele, autres que les amplitudes complexes et la puissance 
de bruit. 

15 Le procede selon I'invention s'applique notamment dans un 

systeme de telecommunications UMTS comportant un reseau d'antennes. 
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Input Multiple Output). 

Dans le cas d'un traitement radar, ces techniques s'appliquent egalement. 
Elles peuvent etre utilisees dans un contexte multi-cible ou mono-cible avec 
echo. 

5 De meme dans le cadre de la radio-altimetrie, plusieurs 6chos du 

signal emis peuvent etre recus et tous necessitent d'etre estimes. 

Sans sortir du cadre de I'invention, le proc^de peut tenir compte 
de I'effet Doppler. Pour cela, une maniere consiste par exemple a estimer les 
10 vitesses en minimisant une fonction multidimensionnelle integrant tous les 
parametres du modele, autres que les amplitudes complexes et la puissance 
de bruit. 

Le procede selon I'invention s'applique notamment dans un 
15 systeme de telecommunications UMTS comportant un reseau d'antennes. 
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REVENDICATIONS 



1 - Precede pour estimer un ou plusieurs parametres d'un canal de 
propagation avec une connaissance a priori du signal dans un systeme 
comportant un ou plusieurs capteurs caractdrise en ce qu'il comporte au 
. moins les etapes suivantes : 
• correler le signal regu par ie ou les capteurs s(t) avec un signal connu . 



• selectionner un nombre d'echantillons du signal obtenu apres I'etape de 
correlation, 

• determiner les parametres du canal tel que x et/ou 9 qui permettent de 
reconstruire au mieux les signaux regus en utilisant une methode au 
maximum de vraisemblance. 

2 - Procede selon la revendication 1 caracterise en ce que les 
caracteristiques du systeme de capteurs sont connues et en ce qu'il 
comporte 

• une etape de correlation avec un signal connu c(t) egal a 1 , 

• les signaux regus sur I'antenne s'exprime sous la forme X=S(0/t)h+B 

• les estimees des parametres x et 0 s'expriment sous la forme 



ou n x s est le projecteur orthogonal a Timage engendree par les vecteurs 
colonnes de S(0,t) 



c(t), 



0 y 9 = argmin||n^(^r)X|f 
= argminlX^^.i-jX} 
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REVENDICATIONS 

5 1 - Precede pour estimer un ou plusieurs parametres cTun canal de 
propagation avec une connaissance a priori du signal dans un systeme 
comportant un ou plusieurs capteurs caracterise en ce qu'il comporte au 
moins les etapes suivantes : 

• correler le signal re?u par le ou les capteurs s(t) avec un signal connu 
10 c(t), 

• selectionner un nombre d'echantillons du signal obtenu apres I'etape de 
correlation, 

• determiner les parametres du canal tel que t et/ou 9 qui permettent de 
reconstruire au mieux les signaux regus en utilisant une methode au 

15 maximum de vraisemblance. 

i 

2 — Procede selon la revendication 1 caracterise en ce que les 
caracteristiques du systeme de capteurs sont connues et en ce qu'il 
comporte - 

20 • une etape de correlation avec un signal connu c(t) egal a 1 , 

• les signaux re?us sur Tantenne s'exprime sous la forme X=S(9,T)h+B " 

• les estimees des parametres t et 9 s'expriment sous la forme 

0,r = argmin||n^(^r)X|| 2 
25 = argminfXtn^fl.rJX} 

ou n x s est le projecteur orthogonal a Timage engendree par les vecteurs 
colonnes de S(9,t) 
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{les caracteristiques du systeme de capteurs sont par example la reponse a 
I'antenne) 

3 - Procede selon la revendication caracterise en ce qu'il comporte une 
5 etape de determination des amplitudes complexes h de la reponse 
impulsionnelle du canal de propagation a partir des estimees des parametres 
estimesx etG. 

4 -Procede selon la revendication 1 caracterise en ce que les 
10 caracteristiques du systeme de capteurs ne sont pas connues et en qu'il 
comporte : - 

• une etape de correlation des signaux regus par le reseau de capteurs 
avec un signal connu c(t) egal a 1 , 

• une etape oli les signaux regus s'expriment sous une forme concatenee 
15 Y ^(tJoc+N ou *F(t) est egal au produit convolue de la matrice unite l N 

avec la matrice S{t) = [s l (rj);- • : .s 1 ^), • • • ,s(r^c;)i et acontient les 
reponses des trajets des differents utilisateurs, 

• une etape d'estimation du vecteur retard x a partir de 



ou II 1 * est le projecteur orthogonal a Timage engendree par les vecteurs 

* 

lignes de 4*(x). 

25 

5 - Procede selon la revendication 1 caracterise en ce qu'il comporte une 



etape de correlation des signaux avec un signal c(t) different de 1, en ce que 



20 



9 = axgmin||ni(T)Y|| 
= argmintr(Y t ni(T)Y) 

T 



les caracteristiques du systeme de capteurs sont connues et en ce qu'il 
comporte une etape d'estimation des parametres x et 9 a partir de 



30 
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(les caracteristiques du systeme de capteurs sont par exemple la reponse a 
Tantenne) 

3 - Procede selon Tune des revendications 1 ou 2 . caracterise en ce qu'il 
comporte une etape de determination des amplitudes complexes h de la 
5 reponse impulsionnelle du canal de propagation a partir des estimees des 
parametres estimes x et 9. 

4 -Procede selon la revendication 1 caracterise en ce que les 
caracteristiques du systeme de capteurs ne sont pas connues et en qu'il 
comporte : 

10 • une etape de correlation des signaux re?us par le reseau de capteurs 
avec un signal connu c(t) egal a 1, 

• une etape ou les signaux regus s'expriment sous une forme concatenee 
Y =4 / (t)cx+N ou m{x) est egal au produit convolue de la matrice unite\l N 
avec la matrice S(r) = [s 1 (ti),--: ,s 1 (rp 1 ),-- ,s(r^«/)] et acontient les 

15 reponses des trajets des differents utilisateurs, 

• une etape d'estimation du vecteur retard x a partir de 

9 = argmm||ni(r)Y|| 2 
= axgniinir (Y f ni(r)Y) 

ou est le projecteur orthogonal a Timage engendree par les vecteurs 
20 lignes de ^(x). 

5 - Procede selon la revendication 1 caracterise en ce qu'il comporte une 
etape de correlation des signaux avec un signal c(t) different de 1, en ce que 
les caracteristiques du systeme de capteurs sont connues et en ce qu'il 
comporte une etape d'estimation des parametres x et 9 a partir de 
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0 T = axgminX w R^ 1 ni(e,r)X w 
ni = I-*(«,r)(» t («,T)Ri i *(»,r))" 1 * t (®.'*-)Ri 1 . 

6 - Procede selon la revendication 1 caracterise en ce qu'il comporte une 
etape de correlation des signaux recus avec c(t) different de 1 et en ce que 
les caracteristiques du systemede capteurs etant inconnues, ^estimation du 
vecteur retard s'exprime a I'aide de 

■ • - 

t = argminY w R; 1 Il£(T)Y w 

; n£ = i - s(r) (swr;^))" 1 smr; 1 

7 - Procede selon Tune des revendications 1 a 6 caracterise en ce qu'il est 
applique dans des applications de type MIMO (Multiple Input Multiple output) 
ou SIMO (Single Input Single Output). 

8 - Dispositif pour estimer un ou plusieurs parametres d'un canal de 
propagation avec une connaissance a priori du signal dans un systeme 
comportant un ou plusieurs capteurs caracterise en ce qu'il comporte au 
moins : 

• un dispositif adapts & corr§ler le signal regu par le ou les capteurs s(t) 
avec un signal connu c(t), 

• un dispositif adapte a selectionner un nombre d'echantillons du signal 
obtenu apres Petape de correlation, et 

• un dispositif adapte a determiner les parametres du canal par une 
methode au maximum de vraisemblance. 

9 - Recepteur de radiocommunication comportant les caracteristiques du 
dispositif selon la revendication 8. 
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0,x = argminX w Rb 1 ni(0,T)X w 
ni = I - t) (&(9, t)R^^(9, t)) _I r)R b -' 

« 
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6 - Procede selon la revendication 1 caracterise en ce qu'il comporte une 
etape de correlation des signaux recus avec c(t) different de 1 et en ce que 
les caracteristiques du systeme de capteurs etant inconnues, Pestimation du 
vecteur retard s'exprime a I'aide de 

10 

t = argminY u, R- 1 n^(r)Y t ' 

T 

n^- = I - §(t) (SC-rjR-^C-r))- 1 S^R" 1 
15 7 - Procede selon I'une des revendications 1 a 6 caracterise en ce qu'il est 

* 

applique dans des applications de type MIMO (Multiple Input Multiple output) 
ou SIMO (Single Input Single Output). 

8 - Dispositif pour estimer un ou plusieurs parametres d'un canal de 
20 propagation avec une connaissance a priori du signal dans un systeme 
comportant un ou plusieurs capteurs caracterise en ce qu'il comporte au 
motns : 

• un dispositif adapte a correler le signal recu par le ou les capteurs s(t) 
avec un signal connu c(t), 

25 • un dispositif adapte a selectionner un nombre d'echantillons du signal 
obtenu apres Petape de correlation, et 

• un dispositif adapte a determiner les parametres du canal par une 
methode au maximum de vraisemblance. 

30 9 - Recepteur de radiocommunication comportant les caracteristiques du 
dispositif selon la revendication 8. 
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